A finite-difference algorithm for full waveform teleseismic tomography

S U M M A R Y We adapt a 2-D spectral domain finite difference waveform tomography algorithm previously used in active source seismological imaging to the case of a plane wave propagating through a 2.5-D viscoelastic medium in order to recover P and S wave speed variations from body waves recorded at teleseismic distances. A transferable efficacy that permits recovery of arbitrarily heterogenous models on moderately sized computers provides the primary motivation for choosing this algorithm. Synthetic waveforms can be generated either by specifying an analytic solution for a background plane wave in a 1-D model and solving for the source distribution that would produce it, or by solving for a scattered field excited by a plane wave source and then adding the background wavefield to it. Because the former approach typically involves a concentration of sources at the free surface, the latter tends to be more stable numerically. We adapt a gradient approach to solve the inverse problem to maintain tractability; calculating the gradient does not require much more computational effort than does the forward problem. The waveform tomography algorithm can be applied in a straightforward way to perform receiver function migration and traveltime inversion.